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The RXTE observations of Scorpius X-l in 1996-1999 are presented. The properties of its quasi-periodic X- 
ray oscillations are studied in detail. The results obtained are used for analysis in terms of the transition-layer 
model (TLM) and the relativistic-precession model (RPM) for a slowly rotating neutron star. Theoretical 
predictions of the two models are compared and their self-consistency is verified. The tilt of the magneto- 
sphere to the accretion-disk plane, the neutron-star mass, and its angular momentum are determined in the 
model frameworks. 
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1. Introduction 

Scorpius X-l is a low-mass X-ray binary with an accreting neutron star and one of the brightest 
X-ray sources. Its spectral properties suggest that Scorpius X-l belongs to the class of Z-type 
sources (Hasinger and van der Klis 1989), whose characteristic feature is a Z-shaped track in the 
color-color diagram. In this interpretation, the spectral properties are presented in the hard-soft 
color indices, each of which is defined as the harder-to-softer flux ratio in the corresponding energy 
band. The Z-shaped track is commonly divided into three parts called branches: the horizontal 
(HB, the upper part of the diagram), normal (NB, the intermediate part), and flaring (FB, the lower 
part) branches. The position on the Z track is generally believed to be associated with the rate of 
accretion in the direction from HB to FB. Six sources are currently known to exhibit Z tracks in the 
color-color diagram: Scorpius X-l, Cygnus X-2, GX 17+2, GX 5-1, GX 340+0, and GX 349+2. 

The power-density spectra (Fourier transforms of the flux) of Z-type sources exhibit low-frequency 
(5-100 Hz) X-ray quasi-periodic oscillation (QPO) peaks. The names of the QPOs correspond to 
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the branch with which their origin is identified: horizontal- (HBO), normal- (NBO), and flaring- 
branch (FBO) oscillations. HBOs (15-100 Hz) can also be detected in the NB spectral state. 
However, as one recedes from HB, the statistical significance of the QPO peaks decreases and they 
become undetectable. When moving along the Z track (from HB to FB) in its NB-FB segment, 
a QPO peak detectable in the range 5-20 HZ (NBO/FBO) emerges in the power-density spectra. 
For all the currently known Z-type sources, QPOs were also found in the range 200-1100 kHz (van 
der Klis 2000). Two kHz QPO peaks (u\ and U2 are the lower and upper peaks, respectively) can 
be simultaneously observed with a frequency difference of ~ 200-400 Hz. For a series of successive 
observations (i.e., on short time scales from several hours to several days), an increase in the flux is 
accompanied by an increase in the frequencies of the two peaks, but their difference (Az/ = v<i — v\) 
is not conserved (van der Klis et al. 1997). 

The observations of Scorpius X-l revealed all the QPO types characteristic of the low-frequency (< 
100 Hz) range in the power-density spectrum. At high frequencies, two kHz QPO peaks {y\andv<i) 
were found in the RXTE/PCA data [first detected by van der Klis et al. (1996)]. Here, we study the 
QPO properties in detail. Our results are used for analysis in terms of the transition-layer model 
(TLM) and the relativistic-precession model (RPM). We compare theoretical predictions of the two 
models and verify their self-consistency. The X-ray flux variability of the source at low frequencies 
(0.1-128 Hz) is investigated by taking into account the power-law behavior of the power-density 
spectrum at frequencies above and below the break frequency. 

2. DATA AND OBSERVATIONS 

For our time analysis, we used data from the PCA (Proportional Counter Array) instrument (Ja- 
hoda et al. 1996) onboard the RXTE observatory (Bradt et al. 1993) retrieved from the Goddard 
Space Flight Center Electronic Archive. 

The source Scorpius X-l was observed by the RXTE observatory during eleven series of pointing 
observations (10056, 10057, 10059, 10061, 20053, 20426, 30035, 30036, 30406, 40020, 40706): in 
February and May 1996, in March, April, and August 1997, in January, February, and from May 
until July 1998, and in January and July 1999. The observations of Scorpius X-l over this period 
correspond to three different observational epochs of RXTE/PCA (1, 3, and 4 in the adopted 
classification), for which the boundaries of the PCA energy channels were changed. 

To construct the power-density spectra, we used observational data with a resolution of ~122 or 
~ 244 /Ltm (2 -13 , 2~ 12 s) from the zeroth to 87th or 249th PCA energy channels. These ranges (0-87 
and 0-249) correspond to a detectable photon flux ~ 1.5 - 23.5, ~ 1.9 - 32.3, and ~ 2.0 - 38 keV 
for epochs 1, 3, and 4, respectively, or extend to ~ 60 keV if the entire accessible PCA energy range 
(channels 0-249) was used. 

Of all the observations, we used only those during which the angle between the source direction 
and the Earth's horizon was more than 10°. During the observations of Scorpius X-l, all five 
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proportional counters were not always switched on to record events. If the operating condition of 
one of the counters changed during a continuous observation (whose duration did not exceed the 
duration of one orbit and was, on average, 3 — 3.5 x 10 3 s), then the time interval during which the 
total count rate changed abruptly was excluded from our analysis. 

To analyze the variability of Scorpius X-l, we constructed its power-density spectra (van der 
Klis 1989) in the frequency range 0.03125-2048 Hz (see Fig. 1). The properties of the low- and high- 
frequency flux variability were investigated in the bands 0.1-256 and 256-2048 Hz, respectively. 
No correction was made for the background radiation and dead time. 
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Fig. 1. The power-density spectrum for the 
Z-type source Scorpius X-l. The data were 
averaged over four consecutive observing in- 
tervals from 05 : 04 : 13 UTC on May 25, 
1996. The total observing time is ~ 10 4 s [see 
a similar power-density spectrum in Titarchuk 
et al. (1999)]. The X-ray QPO peaks are 
shown: HBOs (vhbo) and their second har- 
monic (v^hbo) and the kHz QPO peaks (y\ and 
V2 ) • The break frequency is denoted by VBreak ■ 
A model fit to the data (for more details, see the 
text) is indicated by the dotted line. Broadband 
noise was detected at a statistically significant 
level in the average power-density spectrum at 
low frequencies (vlf)- 



We analyzed individual observations with a duration up to ~ 3.5 x 10 3 s. Fitting the power-density 
spectra by a constant and by a power law at frequencies below and above the break frequency did 
not yield acceptable results (according to the \ 2 test). The main reason was the absence of a sharp 
break and the resulting uncertainty in its measurement. The model in which at frequencies much 
higher [yjvbreak 3> 1) and much lower (v/vbreak *C 1) than the break, each part of the spectrum 
could be fitted by its own power law and the transition between them was not jumplike (previously, 
this model was successfully used for a time analysis of the Z-type source Cygnus X-2; see, e.g., 
Kuznetsov 2001), proved to be more suitable: 

P(v)=A»- a ll + (»/» hieak f]-\ (1) 

The power-density spectra (see Fig. 1) were fitted by this model in the 0.1-128-Hz band with the 
additional introduction of one or two Lorenz lines to allow for the QPO peaks and their harmonics. 
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To take into account the PC A dead-time effect, which causes the overall level to be shifted to the 
negative region [because of this effect, the Poissonian noise level subtracted from all spectra differed 
from 2.0, in units of the Leahy normalization; see van der Klis (1989) and Vikhlinin et al. (1994) 
for more details], we added a constant to the general model. 

In searching for the kHz (in the range ~ 500 — 1200 Hz) QPO peaks v\ and ^2, we analyzed the 
power-density spectra at high frequencies. A constant with the addition of one or two Lorenz lines 
was used as the model. Individual observations in which the detection confidence level of kHz QPO 
peaks exceeded 3d were used for the subsequent analysis. 

Our analysis of the low-frequency flux variability in Scorpius X-l revealed QPOs at a high confidence 
level in all the individual observations that corresponded to the HB or NB spectral states. Note 
that the HBO frequency is one of the main frequencies in most models, including the TLM and 
RPM. The two models (see below) unequivocally establish the correspondence between the three 
QPO peaks: the HBOs, lower, and upper kHz QPOs. 

To obtain the most accurate model parameters (the angle 5 for the TLM and the neutron-star 
mass Mns and angular momentum a for the RPM), which are invariants in each of the theories 
under consideration, it was necessary to calculate the three simultaneously observed main QPO 
frequencies (i.e., vhbo-, v \-> v 2\ see Fig. 1) as accurately as possible. Therefore, of all the preselected 
data in which both kHz QPO peaks (y\ and v<i) were found, we used only those in which the main 
HBO peak (vhbo) and its second harmonic (U2HB0) were detected at a confidence level above 4a. 
In such observations, the two kHz QPO peaks were also reliably detected. 

3. LOW-FREQUENCY X-RAY FLUX VARIABILITY 

3.1. The Power of the HBO peaks 

In Fig. 2, rms is plotted against HBO frequency. The flux variability (i.e., rms 2 ) at the peak 
corresponds to the integral of the Lorenz line that fits best the QPO peak in the power-density 
spectrum. For the Z-type sources, to which Scorpius X-l belongs, an increase in HBO frequency 
is accompanied by a decrease in rms (see, e.g., van der Klis 2000). As a result, the detection 
confidence level of HBO peaks is lower at high frequencies (~ 50 — 60 Hz). Figure 2 shows the 
range of HBO frequencies found for Scorpius X-l. During all observations, the deviation from the 
mean vjjbo ~ 43.5 Hz did not exceed ~ 5 Hz (i.e., it changed by no more than 10% during all 
observations) . 
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Fig. 2. The rms of the QPO peaks versus 
HBO frequency. Individual observations are 
presented. As rms decreases, the detection con- 
fidence level of the QPO peaks falls, reaching 
<~ 4cr for rms in the range <~ 0.4 — 0.5%. 



3.2. Horizontal- Branch Oscillations and Their Harmonics 

In Scorpius X-l, quasi-periodic HBO peaks were found at low frequencies (< 100 Hz). Together 
with the oscillations at the main frequency vhbo-, we also found QPO peaks at frequencies close to 
the multiple ones. Note that the detection confidence level of HBO peaks decreases with frequency 
(see Fig. 2). That is why no multiple peaks were found (at a confidence level > 4<r) for which 
vhbo > 48 Hz (the corresponding oscillation frequency of the second harmonic is v 2 HBO > 96 Hz). 

The frequency ratio r 2 i\ = ^hbo/^hbo of the multiple HBO harmonics is shown in Fig. 3. We 
clearly see that r 2 i\ differs from 2.0: only a third of the 39 data points are above 2.0 (the dotted 
line in Fig. 3). Fitting by a constant yields r 2 /i = 1.965 ± 0.004, Xred = 3-45 (in what follows, 
Xred = X 2 /d.o.f is the reduced \ 2 value per degree of freedom), while fitting by a constant with 
a break yields r 2 /i = 1.976 ± 0.004, \l e d = !- 05 ( anci r 2/i = -0.05 x (u HB o ~ 45.31) in the range 
of a linear dependence). When the scatter of data points exceeds the stochastic scatter and when 
the Xred vame is much larger than unity (as in the case of a constant fit), it should be assumed 
that the uncertainty could be produced by stochastic processes rather than by the statistics alone. 
We can then ignore the measurement error and determine the mean f 2 /i = jj J2 r 2/i an d the rms 

deviation from it a = "sj ^^'^Ji^ 11 ^ ■ I n this case, we obtained the following ratio of the second 
and first harmonics: f 2 /i = 1.973, a = 0.050. Under such an assumption, the harmonic ratio is 
compatible with 2.0, within the error limits. 



-6- 



Fig. 3. The ratio of multiple HBO harmon- 
ics versus frequency vhbo- The solid, dash 
dotted, dashed, and dotted lines represent, re- 
spectively, a constant fit to the data with sta- 
tistical errors (minimization according to the 
X 2 test), a constant fit with a break, the mean 
of individual observations (without errors), and 
the assumed harmonic ratio 2.0. 
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Let us consider in more detail the case where the scatter of r 2 /i values about the mean is deter- 
mined by the statistics alone and is incompatible with 2.0, within the error limits. Resonances in 
nonlinear oscillatory systems can emerge at frequencies at which an external periodic force acts on 
the oscillator. This requires that its frequency 7 satisfy the condition 

7 = ^0-, (2) 

where p and q are integers and ujq is the cyclic frequency of free system oscillations. 

The power of the main peak must be maximal among all of the harmonics, because the intensity 
of resonance phenomena rapidly decreases with increasing p and q (Landau and Lifshitz 1965). In 
addition, if the frequency of the periodic external force 7 differs from ujq by a small value e (in the 
simplest case, p = 1 and q = 1), then in nonlinear oscillations, this can lead to a displacement of the 
rms maximum from the oscillator resonance frequency (Landau and Lifshitz 1965). In Scorpius X-l, 
the largest rms among all the HBO harmonics corresponds to the frequency vhbo ( see Fig- 2). We 
assume from the outset that the latter (and this assumption proved to be valid) was the fundamental 
harmonic and the deviation of the frequency ratio of the second and first harmonics from 2.0 is 
admitted in the forced oscillations of a nonlinear oscillator under the action of a periodic force. 

Note that in this case, the frequency of the observed low- frequency oscillation peak is vrbo 7^ 
cjo/27t. This is particularly important in comparing the self-consistency of the theoretical models 
under consideration, which interpret the emergence of QPOs (to be more precise, the following 
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three QPOs: vhbo-, u i, and ^2) in binary systems with neutron stars. 

If the deviation of the ratio of the multiple HBO harmonics v 2 /i from 2.0 is assumed to be at- 
tributable to forced nonlinear oscillations, then the observed vhbo must be close but not exactly 
equal to ujq/2tt, the assumed frequencies of free oscillations of the oscillator. Note that lvo/2it rather 
than vhbo corresponds to the low-frequency QPOs in the models considered below. 

For this reason, for our analysis in terms of each of the models considered below, we use two values 
for the low- frequency QPOs: vhbo an d vhbo- These values are the fundamental HBO harmonic 
and the mean between half the frequency of the second harmonic v 2H bo and the fundamental har- 
monic vhbo- F° r the latter case, from all the observations, we chose those in which the confidence 
level of each HBO peak was higher than 4cr. 

4. THEORETICAL MODELS 

For our analysis in terms of the two models in question, we used the average spectra and the main 
frequencies: vhbo (as well as vhbo = ^/2(vhbo + l/2^2_ffSo)) the mean between the fundamental 
harmonic and half the second harmonic) , v\ , v 2 ■ The following approximate asymptotic relations 
between the frequencies are known from the observations of low-mass X-ray binaries (Stella et al. 
1999): (1) Av = v 2 - v x oc z/f ; and (2) vhbo °t v \- These observed features are consistent with 
each of the theories under consideration. 

4.1. The Transition-Layer Model (TLM) 

We consider the motion of a clump of matter, the QPO source on the accretion-disk surface, in 
a Keplerian orbit around a neutron star (NS) in terms of this model. The magnetospheric axis 
is assumed to be not aligned with the normal to the disk surface but makes an angle 5. After 
multiple passages through a slightly tilted magnetosphere, the clump comes under the effect of 
Coriolis forces. These forces cause the main Keplerian oscillation frequency vk to split up into two 
oscillation modes: radial (v^) and perpendicular to the disk plane (vl). The two modes are the 
solution of the equation for the rotation of a body in a noninertial frame of reference [see Osherovich 
and Titarchuk (1999) for more details on its derivation and solution]. In this model, the lower and 
upper kHz QPO peaks correspond to the Keplerian (vk = v\) and hybrid (v^ = v 2 ) frequencies, 
respectively. The relation between the kHz QPO peaks is given by 

v^ivl + tn/Tr) 2 ] 1 / 2 , (3) 

where f2 is the angular velocity of the magnetosphere. The oscillation mode perpendicular to the 
disk surface is defined as 

v L = (M/tt) (v K /v h )sin 8. (4) 
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Since all three characteristic frequencies vl, "k, and Vh are known from observations, the sought-for 
angle 5 can be calculated by using (4): 

5 = arcsin[(i^ - i>K)~ ll2 {v L v h /v K )\. (5) 

To a first approximation, the angular velocity of the magnetosphere Q is constant. A more accurate 
equation that describes the dependence of O on radius can be derived in the multipole magnetic- 
field approximation (see Osherovich et al. 1984). Assuming the contribution of the quadrupolc 
component to the magnetic-field strength in the equatorial plane to be negligible and taking into 
account the dipole and octupole components, the final equation for the angular velocity can be 
represented as (Osherivich and Titarchuk 1999) 

n/2ir = C + C\u 4 J 3 + c 2 4 /3 + C 3 z4, (6) 

with C 2 = -2(CiC 3 ) 1/2 - 

To reconstruct the magnetospheric profile for Scorpius X-l, according to Eq. (3), it will suffice to 
have data only for the high-frequency range of the flux variability, i.e., v\ and v-i- Therefore, we 
used all the data in which kHZ QPO peaks were found by removing the additional condition that 
the main HBO peak vhbo and its second harmonic vihbo should be detected in the power-density 
spectrum at low frequencies at a confidence level above 4er. 

The angular velocity is plotted against Keplerian frequency in Fig. 4. Using (6), we managed to 
obtain the magnetospheric profile suggested by the model under consideration with the following 
parameters: C = 354 Hz, C x = -3.54 x 1(T 2 Hz" 1 / 3 , C 2 = 9.99 x 1(T 6 Hz~ 5 / 3 , and C 3 = 
-7.21 x 10~ 10 Hz" 3 . In fitting the data, we assumed C2 to be a free parameter, whereas in the 
model, C2 = — 2(CiC l 3) 1 / 2 =1.01xl0 -5 . The similar values of C2 determined by two independent 
methods are a weighty argument for the approximation used for the magnetospheric profile. 

Using (5), we obtained 5 for each triplet of frequencies: vhbOi v i-> an d v-i- I n Fig- 5 the HBO 
frequency uhbo is plotted against the angle 5 between the normal to the disk surface and the 
magnetospheric axis. Among all three frequencies used in calculating 5, the strongest correlation 
is observed between 5 and vhbo- F° r this reason, the dependence of 5 on low-frequency QPOs is 
analyzed in the model under consideration (below, we show that this is also true for the relativistic 
precession model), because the model invariants must be conserved, irrespective of the measured 
QPO frequency range. 

Table 1 gives data on the constant fit to 5 in two cases: when vl = vhbo (see Fig. 5) and 
vl = vhbo = 1/2(vhbo + 1/2i / 2H_bo) (see Fig. 6), i.e., with the inclusion of the second HBO 
harmonic. The Xred va l ue is given in each case. The mean angle 5 and its rms deviation od are 
also given in Table. 1. 

To reduce the possible stochastic scatter of data points (which probably dominates over the statisti- 
cal scatter), we performed averaging depending on the HBO frequency. The upper panels in Figs. 5 
and 6 show the data obtained in this way. Whereas the data in Fig. 5 reveal a clear correlation 
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Fig. 4. The assumed angular velocity of the 
magnetosphere in the TLM, in units of Q/2ir, 
versus Keplerian frequency [y K — v{). 
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between <5 and vhbo^ the data points in Fig. 6 virtually lie on the same line. If a larger number of 
points in each bin were averaged (in our case, we combined them by five), then 5 (see Fig. 6, where 
5 is plotted against vhbo) would be even more compatible with a constant, according to the \ 2 
test. 

Let us compare the rms deviation between the HBO harmonic ratio r 2 /i determined above and <5 
(see Table 1). It turns out that o> 2/1 « 2.5%, while as ~ 2.0%. Since, according to (5), 5 oc vhbo 
at small angles and since vubo increasingly deviates from uj/2tt to larger values with decreasing 
r2/i, the correlation for 6 (clearly seen only in Fig. 5) can result from the emergence of nonlinear 
oscillations considered in the previous sections. 



4-2. The Relativistic Precession Model (RPM) 

The high-velocity motion of matter in a strong gravitational field can generate oscillations at- 
tributable to general-relativity effects. The RPM (see, e.g., Morsink et al. 1999) considers the 
motion of a point mass around a gravitating center. If a particle moves in an orbit that does not lie 
exactly in the equatorial plane of a compact object but is inclined at an infinitesimal angle, then 
the particle orbit will precess. 

For a nonrotating gravitating center (i.e., in the Schwarzschild approximation), the expression for 
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Fig. 5. The angle between the magnetospheric 
axis and the normal to the disk surface versus 
HBO frequency vhbo = v l- A constant fit is 
indicated by the solid line. The lower panel 
shows individual observations of Scorpius X-l 
with a duration up to ~ 3.5 ks. The upper 
panel shows the same observations as in the 
lower panel averaged over five points according 
to the frequency vhbo- 
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the particle angular velocity matches the classical formula for Keplerian motion. However, for a 
NS or a black hole with an intrinsic angular momentum, the azimuthal frequency can be expressed 
(Bardeen et al. 1972; Stella et al. 1999) in a system of units with G = 1 and c = 1 as 



2vr 1 + a 




-l 



(7) 



Below, we consider only the case where the compact object and the test particle corotate: v§ > 
0, a > (a is the relative angular momentum). The epicyclic frequency v r together with the 
azimuthal frequency determine the orbital periastron rotation, v per = v^ — Vr, while the frequency 
in a direction perpendicular to the disk plane vq determines the nodal precession, v no d = \ v 4>~ v e\- 
The corresponding equations were derived by Okazaki et al. (1987) and Kato (1990): 



u r = 



M IM a* 

6 1- 8a\l — - 



vq = 



M a 2 
4cm/ — ^ + 3-^-. 



(8) 



(9) 



In contrast to the TLM, the Keplerian rotation frequency in the RPM corresponds to the upper 
QPO peak, i>2 = u^, while the observed periastron precession frequency corresponds to the lower 
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Tabl. 1. The angle between the normal to the 
disk surface and the magnetospheric axis for var- 
ious HBO values. 



Parameter 


VL = VHBO 


V L = v a HBO 


5 


5.612 ±0.007 


5.545 ± 0.006 


•^red 


5.53 


2.94 


5 b 


5.600 ±0.114 


5.555 ± 0.086 


N, bin 


39 


39 



a VHBO = l/2(v H BO + l/2ls 2H Bo)- 

b The quantity S is the mean of all S with the error <td 
equal to the standard deviation. 



QPO peak, v\ = v per = v§ — u r . In the Kerr approximation, ^ fg (for the Schwarzschild case 
with a = 0, the frequencies are equal: = isg), the HBOs are identified with nodal precession: 
v no d = 1 1/0 — vq\. Some authors (e.g., Stella et al. 1999) believe that vhbo is an even harmonic of 
f riod- In contrast to vhbo = ^nod, this relation causes the mass of the compact source in Eqs. (7)-(9) 
required by the RPM to decrease. 

In the RPM approximation under consideration, the NS mass and relative angular momentum are 
invariants and do not depend on Keplerian frequency. Thus, we can determine Mjys and a from 
Eqs. (7)-(9) using the three observed main frequencies: vhbo-, v \-> an d V2- 

As in the section devoted to the TLM, here, the invariant (more specifically, the NS mass M^s 
suggested by the model) was also determined in two cases: v n od = ^hbo an d v n od = ^hbo = 
\j1{vHBO + l/2f2ffso) (he., with the inclusion of the second HBO harmonic). Since the derived 
M^s proved to be most dependent (as in the TLM) on low-frequency HBOs, Figs. 7 and 8 show 
Mns as a function of vhbo- 



Tabl. 2. The NS mass and its relative angular momentum for various HBO values. 



Parameter 


Vnod = VHBO a 


^nod = VHBO 


2v no d = vhbo 


2v no d = UHBO a 


M NS ,M Q 


2.680 ± 0.002 


2.675 ± 0.002 


2.339 ± 0.002 


2.338 ± 0.002 


Xred 


14.3 


16.3 


16.4 


17.5 


M b NS ,M Q 


2.677 ±0.050 


2.674 ±0.053 


2.341 ± 0.045 


2.341 ± 0.045 


a c , 


1367 ± 2 


1358 ± 2 


663 ± 1 


658 ± 1 


Xred 


8.2 


13.7 


8.2 


13.6 


N, bin 


39 


39 


39 


39 



a V H BO = l/2(v H BO + l/2u 2 HB0)- 

b M N $ corresponds to the mean of all M N s with the error equal to the standard deviation 
<Jd- 



c a = J/(c x Mjvs). 

Table 2 gives a constant fit to the data (shown in Figs. 7 and. 8). We clearly see a large Xred va lue 
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Fig. 6. The angle between the magnetospheric 
axis and the normal to the disk surface ver- 
sus HBO frequency vhbo = v l — 1/2(vhbo + 
1/2i>2hbo)- A constant fit is indicated by the 
solid line. The lower panel shows individual ob- 
servations of Scorpius X-l with a duration up 
to ~ 3.5 ks. The upper panel shows the same 
observations as in the lower panel averaged over 
five points according to the frequency vhbo- 
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for each 

Vnod- compared to Table. 1, which gives a fit to 5 (in the TLM), the corresponding Xred 
values are larger by a factor of ~ 3 — 6. Nevertheless, the rms deviation ctm ns ls a l so ~ 2.0%, as is as 
in the TLM. In contrast to Fig. 4, in which S is plotted against vhbo and is essentially compatible 
with a constant, we clearly see a linear correlation of M^s with vhbo 111 Fig- 8: as vhbo increases 
from ~ 39.5 Hz to ~ 46.5 Hz, the required M^s decreases from ~ 2.8M© to ~ 2.6M . 

The relativistic precession model makes it possible to estimate the compact-source mass and to 
determine the Keplerian orbital radius. For comparison, Fig. 9 shows the mass-radius relations for 
various equations of state for neutron stars taken from Miller et al. (1998). 

4-3. Comparison between the Assumed Invariants of the Models under Consideration 

To check whether the two models are self-consistent and to determine whether the derived invariants 
5 and M^s correlate with the HBO frequency (see Fig. 5-9), we used Spearman's nonparametric 
correlation test. The coefficient r s equal to unity in magnitude points to a total correlation between 
the two quantities. In each case under consideration, we give the probability (prob) that the 
correlation found is the result of statistical fluctuations. This test is also convenient in the case 
where it is necessary to establish which of the several correlations under consideration is strongest. 
Our estimates are given in Tables 3 and 4. Note that the data in Table 4 correspond to the case 
where the statistical scatter was slightly reduced through averaging. Thus, the RPM model (see 
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Fig. 7. NS mass versus observed HBO fre- 
quency i^h bo- The solid line indicates a con- 
stant fit to the data. The lower panel shows 
individual observations of Scorpius X-l with a 
duration up to ~ 3.5 ks. The upper panel shows 
the same observations as in the lower panel 
averaged over five points according to the fre- 
quency V H BO- 
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Table 4) is inconsistent (r s = —1; the total negative correlation between M^s and vhbo) and the 
assumed NS mass is not an invariant in this model. 

Allowance for the second harmonic in determining the low-frequency QPO peak (more specifically, 
vhbo = 1/2(vhbo + l/2^2HBo)) results in a reduction of the correlation coefficient r s for 5 in the 
TLM, while the reverse is true for Mns- In particular, we see from Fig. 8 that the data points are 
well fitted by a straight line: x 2 re d = 4-31 (see Table 2, where Xred = 16.3 for the same case but 
for a constant fit to the data). Figures 6 and 8 (upper panels) show that averaging can reduce the 
stochastic scatter when determining the model invariants. However, for the TLM (see Fig. 6), S 
is virtually constant, 5 = 5.55, with the rms deviation a$ = 0.03 (which is ~ 0.5%), while for the 
RPM (see Fig. 8), Mjvs strongly correlates with vhbo an d is incompatible with a constant (the 
rms deviation is ^m ns ~ 1.7%, which roughly corresponds to om ns \ see Table 2). 

5. DISCUSSION AND RESULTS 

Both theoretical models suggest the existence of invariants. These invariants are the tilt of the 
magnetosphere to the disk normal S for the TLM and the NS mass M^s and its relative angu- 
lar momentum a for the RPM in the Kerr approximation. Since the RPM did not consider the 
quadrupole moment of the compact source and disregarded its oblateness (which can give a con- 
tribution of ~ 10 — 15% to the nodal precession frequency), we managed to avoid the inclusion of 
various theoretical equations of state for neutron stars in our analysis. However, it was of interest to 
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Fig. 8. NS mass versus observed HBO 
frequency v nod = T> m HBO = l/2(v HB o + 
l/2v2HBo)- The solid line indicates a constant 
fit to the data. The lower panel shows individ- 
ual observations of Scorpius X-l with a dura- 
tion up to ~ 3.5 ks. The upper panel shows the 
same observations as in the lower panel aver- 
aged over five points according to the frequency 
vhbo- 



compare our results with the constraints imposed on the equations of state. Figure 9 shows some 
of these constraints taken from Miller et al. (1998). The curves correspond to the mass-radius 
relation for nonrotating neutron stars. Although the angular velocity for Scorpius X-l is nonzero, 
the correction for spin (here, v s « 200 — 400 Hz) can be disregarded. 

The assumed NS mass Mjys and the Keplerian orbital radius R of the particle rotating around it (see 
Fig. 9) derived from Eqs. (7)-(9) must be consistent with the constraints imposed by the equations 
of state. More specifically, M^s < Meos, R > Reos (here, Meos and Reos are the NS mass and 
radius for various equations of state). Only the equation of state "L" suggest that a massive NS (in 
our case, ~ 2.7M Q ) could exist and it agrees with the derived m the two cases [y n od = v hbo 
and 2v no d = vhbo)- At the same time, the other two equations of state ("A" and "UU") are 
incompatible with the RPM results in the Kerr (i.e., with NS rotation) approximation. Moreover, 
we obtained M/vs for Cygnus X-2 that was compatible, within the error limits (M^s = 2.69 ±0.10 
for v no & = vhbo and M^s = 2.29 ± 0.09 for 2v no d = vhbo] Kuznetsov 2002), with M^s for 
Scorpius X-l (see Table 2). Such a coincidence is not an argument for the relativistic precession 
model either. 

In addition, we showed that a nonparametric test for correlation between the assumed NS mass in 
the RPM and the HBO frequency yields a positive result. Thus, the RPM is not a self-consistent 
model in which the invariants are conserved and the derived mass of the compact object is large 
enough for the existing equations of state for neutron stars. At the same time, for the transition- 
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Fig. 9. Relationship of the NS mass to the 
nodal precession frequency and Keplerian or- 
bital radius. The circles correspond to the 
mass-radius relation in the case where the 
nodal precession frequency is identically equal 
to the HBO frequency umbo] the diamonds cor- 
respond to the case where it is assumed to be its 
even harmonic (for more details, see the text). 
The solid lines [taken from Miller et al. (1998)] 
represent the mass-radius relations for several 
equations of state for neutron stars. 
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layer model, the derived value is compatible with a constant under certain conditions (e.g., using 
the theory of nonlinear oscillations; for more details, see Kuznetsov and Titarchuk 2002). If the 
uncertainty in the calculated angle 5 does not exceed ~ 2% (which is lower than the rms deviation 
for 5 in Scorpius X-l; see Table 1), then its value may be assumed to be constant. 

The observed inconsistency of the RPM may stem from the fact that the orbital inclination of the 
test particle in the model simplification is assumed to be infinitesimal. An exact solution for an 
arbitrary inclination is given in Sibgatullin (2001). In such an analysis, the frequency uhbo Ior 
fixed mass and angular momentum can change by a factor of ~ 3 as the inclination changes from 
to it/2 (for the marginally stable orbit at M^s = 2M and = 1200 Hz). This dispenses with 
the need to take the nodal precession frequency u no( i to be equal to 1/2 of vhbo to obtain more 
acceptable NS-star masses. 

In addition, it is worth noting that the angle 5 in the RPM may not be strictly an invariant. Higher 
energy release in the disk can affect the curvature of its inner region. The tilt angle 5 between the 
magnetospheric axis and the normal to the disk surface can then differ for different source fluxes. 
As a result, the angles 5 can differ for the same vhbo (see Figs. 5 and 6). 
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Tabl. 3. Correlation between model invariants, S in the TLM and Mms m the 
RPM, and HBO frequency. 



Parameter 


Model 


TLM 


RPM 


VL = V H BO 


VL = V H BO a 


V no d = VHBO 


"nod = l ? HBO a 


r s 


0.567 


0.359 


-0.715 


-0.789 


Prob 


2 x 10" 4 


3 x 10" 2 


3 x 10~ 7 


2 x 10~ 9 


N, bin 


39 


39 


39 


39 



a V H BO = l/2(v H BO + l/2v 2 HBo) 



Note. We used the individual observations shown in the lower panels of Figs. 5-8. 
The Spearman coefficient r s close to f or —1 corresponds to a total positive or 
negative correlation between the two quantities. The probability that the observed 
correlation is the result of purely statistical fluctuations is given for each case. 



Tabl. 4. Correlation between model invariants, 8 in the TLM and Mns in 
the RPM, and HBO frequency. 



Parameter 


Model 


TLM 


RPM 


VL = VUBO 


VL = VHBO a 


"nod = VUBO 


v nod — u HBO 


r s 


0.738 


0.595 


-0.929 


-1.000 


Prob 


3.7 x 10" 2 


1.2 x lO" 1 


8.6 x 10" 4 


0.0 


N, bin 


8 


8 


8 


8 



a V U BO = l/2(v H BO + l/2v 2 HBo) 



Note. We used the data (averaged over four or five observations) shown in the upper 
panels of Figs. 5-8. The Spearman coefficient r s close to f or — f corresponds to a 
total positive or negative correlation between the two quantities. The probability 
that the observed correlation is the result of purely statistical fluctuations is given 
for each case. 
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